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We present a comprehensive numerical study of a microscopic model of the fractional quantum 
Hall system at filling fraction v — 5/2, based on the disc geometry. Our model includes Coulomb 
fH | interaction and a semi-realistic confining potential. We also mix in a three-body interaction in 

I , some cases to help elucidate the physics. We obtain a phase diagram, discuss the conditions under 

which the ground state can be described by the Moore-Read state, and study its competition with 
neighboring stripe phases. We also study quasihole excitations and edge excitations in the Moore- 
Read-like state. From the evolution of the edge spectrum, we obtain the velocities of the charge 
and neutral edge modes, which turn out to be very different. This separation of velocities is a source 
of decoherence for a non-Abelian quasihole/quasiparticle (with charge ±e/4) when propagating at 
the edge; using numbers obtained from a specific set of parameters we estimate the decoherence 
length to be around four microns. This sets an upper bound for the separation of the two point 
contacts in a double point-contact interferometer, designed to detect the non-Abelian nature of 
such quasiparticles. We also find a state that is a potential candidate for the recently proposed anti- 
Pfafhan state. We find the speculated anti-Pfaffian state is favored in weak confinement (smooth 
edge) while the Moore- Read PfafBan state is favored in strong confinement (sharp edge). 
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> 1 I. INTRODUCTION 

in 

Fractional quantum Hall (FQH) liquids are remarkable many-electron systems that possesses nontrivial topological 
i order i Such topological order is reflected in experimental measurable quantities, including the (fractional) charge and 
statistics angle of the gapped quasiparticle/quasihole excitations supported by the system, and the spectra of gapless 
edge excitations. By now a large number of different FQH states, usually labeled by a Landau level (LL) filling factor 
■ v = N e /N<p (where N e is the number of electrons and the number of flux quanta enclosed in the system), have 
been observed experimentally. Most of these FQH states are the Laughlin states or their hierarchy descendents. These 
are Abelian FQH states whose quasiparticles obey Abelian fractional statistics, and their edge excitations are made 
of chiral bosonic modes. The quasiparticle charge has been measured in some of these states^ and measurement of 
the statistics angle has been attempted recently^ Edge excitations of such FQH states have also been probed using 
^ ' electron tunneling^ and other methods. 

Recently much interest and attention have focused on a special FQH state with filling factor v = 5/2, first observed 
twenty years ago£ Interest in this system is driven in part by numerical work^i that suggests the Moore- Read Pfaffian 
stated is likely realized in the half-filled first excited Landau level (ILL) at this filling factor.— The Moore-Read state 
is qualitatively different from the Laughlin states and their descendents, in that it is a non-Abelian FQH state, whose 
quasiparticles obey non-Abelian statistics 10 and whose edge excitations include a branch of fermionic modei^ It has 
been suggested that non-Abelian quasiparticles can be used for topological quantum computation i 12 i 13 ' 14 i 15 ' 16 further 
fueling the interest in this system. It is known that the Pfaffian state is not particle-hole symmetric. It was pointed out 
very recentl y 17 ' 18 that the particle-hole transformed Pfaffian state, termed the anti-Pfaffian state, is also a contender 
at v — 5/2. These two states are closely related, but different in important ways that have experimental consequences. 
While the electron-electron interaction is particle-hole symmetric when projected onto a half-filled LL, which suggests 
these two states would be degenerate if this were the only term present in the microscopic Hamiltonian, in reality the 
degeneracy between the Pfaffian and anti-Pfaffian states is lifted by terms that break particle-hole symmetry. These 
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include LL mixin g i 18 and, as we will show later, confining potential. 

In this paper we report results of a numerical study of the half filled ILL in a disc geometry. Our study is 
complementary to earlier numerical works based on the sphere^ and torus^ geometries, because the disc is the only 
geometry that allows one to study edge states and the closely related physics associated with a confining potential. Our 
results can be briefly summarized as follows. By varying both the electron-electron interaction and confining potential, 
various types of ground states are stabilized. We find that the Moore-Read Pfaffian and possibly the anti-Pffafian 
ground states are realized in different regions of the parameter space of our model. Within our model they appear to 
be separated by an intermediate state that we interpret as a stripe state. We further study the quasihole and edge 
excitations of the Pfaffian state, and show that they indeed have the properties predicted by theory. Furthermore, we 
are able to extract the velocities of the Pfaffian edge modes, which are of importance in addressing both qualitative 
and quantitative issues that arise in experimental studies of the edge states^ especially those involving quasiparticle 
tunneling in a double point-contact interferometeri 20 i 21 i 22 i 23 i 24 i 25 Some of our results were briefly reported in an earlier 
letter ,22 

The rest of the paper is organized in the following way. We describe our microscopic model and its Hamiltonian 
(a mixture of Coulomb interaction and three-body interaction) in Sec. [ill In Sec. IIIII we study various competing 
ground states, which emerge as the lowest energy states in the exact diagonalization study. In Sec. II VI we discuss the 
trapping of charge +e/4 quasiholes by local potentials (generated by, say, an atomic force microscope tip) in certain 
ground states which are supposed to be in the same universality class as the Moore-Read state. We then discuss the 
evolution of the edge spectrum with the variation of interaction in Sec. [V] in particular, we provide an estimate of 
the charge and neutral velocities in a real system based on our model, and discuss the implication in the decoherence 
in double point-contact interference experiments. In Sec. IVII we demonstrate the non-Abelian nature of a charge 
+e/4 quasihole by comparing the edge spectra of a system with and without the quasihole. Potential instability in 
the fermionic edge mode is found. In Sec. IVIII we discuss a potential candidate that emerged from the numerical 
calculations for the recently proposed anti-Pfaffian state, and speculate on its stability conditions. We summarize 
our results in Sec. IVIIII We leave the technical details of the identification of edge states in a system with mixed 
three-body and Coulomb interaction to Appendix [X] The detailed analysis of the evolution of edge states in the pure 
Coulomb limit is presented in Appendix iBl 



II. THE MICROSCOPIC MODEL 



We consider a microscopic model of a two-dimensional electron gas (2DEG) confined to a two-dimensional disk, 
with a mixed Hamiltonian 

H = XH 3B + (1-X)H G . (1) 

Here, the parameter A interpolates smoothly between the limiting cases of a purely three-body Hamiltonian H^b 
(A = 1) and a pure two-body Coulomb Hamiltonian He (A = 0). In the following we measure length in units of the 
magnetic length Is = \/Tic/eB (B is the magnetic field) and energy in units of e 2 /els (e is the dielectric constant), 
such that all quantities that appear later are in units of some combination of the two based on their dimensionality. 
Explicitly, the three-body interaction H^b has the form 

H 3B = - ]T ^, fc [V 2 V 2 (V 2 + V 2 3 )S(r t - rj )S(Ti - r fe )], (2) 

i<j<k 

where S is a symmetrizer: 6123 [/123] = /123 + /231 + /312, where / is symmetric in its first two indices. The iV-electron 
Pfaffian state proposed by Moore and Read 8 - for a half-filled lowest Landau level (OLL) 



^mr(zi,z 2 ,...,z n ) 

i<j K i ) 

is the exact zero-energy ground state of H%b with the smallest total angular momentum Mq = N(2N — 3)/2. In 
Eq. ((3]), the Pfaffian is defined by 

N/2 

PiM v = 2 "l*(N/2)\ £ SgnCT II ^(2fc-iM2fc) (4) 
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for an N x N antisymmetric matrix with elements My. In reality, three-body interaction is present due to finite 
Landau level mixing. The three-body Hamiltonian also has other zero-energy states, known as the edge states, which 
will be discussed in Sec.|V] We note that while the Moore-Read Pfaffian wave function Eq. (|3|) is written for electrons 
in the OLL, it is straightforward to generate the corresponding wave function for electrons in the ILL by applying the 
LL raising operator to every electron. For the rest of the paper we will use the OLL version of various wave functions 
to simplify our discussion, with the understanding that the ILL version of the wave function is generated the same 
way. 

However, there is a more transparent formulation of the three-body Hamiltonian in terms of projection operators^ 
which can be written as 

H 3 b = E E \i'M(hj ) k)){ipM(i,j>k)\, (5) 

M i<j<k 

where ipMih j> k) is a three-particle wave function specified below [Eq. and M is the total angular momentum of 
the state. The Hamiltonian for three particles produces a single non-zero eigenvalue which is unity (provided sufficient 
number of orbitals are allowed) as a true projection operator should. This is the most natural way to define the scale of 
the three-body Hamiltonian. It is simpler to analyze tpM for bosons first. The corresponding expression for fermions, 
as usual, is obtained by multiplication of an appropriate Jastrow factor. The Moore-Read wave function for bosons 
contains one unit of relative angular momentum in the Laughlin factor for each pair, instead of two for fermions. 
As a result when three particles are brought together the relative angular momentum is 2 = 3x1 — 1 (instead of 
5 = 3x2 — 1). We now need to project out all relative angular momenta smaller than two. In this case, the only 
possibility is angular momentum zero (see Ref. [2?] for details). The relative wave function is thus a constant. The 
total angular momentum M will have to be absorbed by the center of mass wave function, which is (zi + z 2 + z 3 ) M . 
For fermions the normalized wave function is: 

->Pm(zi, z 2 , z 3 ) = B M (zi + z 2 + z 3 ) M ~ 3 J(zi,z 2 , 23), (6) 

where J{z\, z 2l z 3 ) = [z\ — z 2 ){zi — z 3 )(z 2 — z 3 ), and the normalization factor is: 



1 / 3M-4 

Bm = (2^)3/2 y 2 M + 2 (Af- 3)!' (7) 

The total angular momentum of the Jastrow factor J is 3 and that of the center of mass is M — 3, giving a total 
angular momentum M. 

The three-body interaction H^b has a rather simple form in the occupation space: 



H ^ = E E ^(M) 



mi>m2>m3 m4<m5<mg 

and 

U({m t }) = V(m 1 ,m 2 ,rn 3 )V(m 4: ,rn 5 ,rn 6 ), (9) 
where V is a completely antisymmetric function of its arguments. With M = mi + m 2 + m 3 we have: 



/ (A/ -1)' 

V(mi,m 2 ,m 3 ) = \ M : : ^Hmifrai - 1)}, (10) 

y 2 x S M mi\m 2 \m 3 \ 

and A is the antisymmetrizer in mi, m 2l and m 3 . The difference between the spectra of this form of the Hamiltonian 
and the one with the (5-functions is an overall factor of 7r 2 /8 in the latter. While there are more efficient ways to 
obtain the Moore-Read stated that avoid diagonalizing a three-body Hamitonian, here we need the H 3 b to generate 
the spectrum of the mixed three-body and the two-body Coulomb Hamiltonian. 

The Coulomb Hamiltonian He includes a two-body Coulomb (1/r) interaction and a one-body confining potential 
provided by the neutralizing background charge distributed uniformly on a parallel disk of radius R, placed at a 
distance d above the 2DEG. This distance parameterizes the strength of the confining potential, which decreases 
with increasing d. The rotationally invariant confining potential comes from the Coulomb attraction between the 
background charge and the electrons. Using the symmetric gauge we can write down the following Hamiltonian for 
the electrons confined to the ILL: 

H G = ^ E V rnn C L+l C i C n+lC m + E U ^ m C mi (11) 
mnl m 
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v = 2 IQHE edge (OLL) 

FIG. 1: (color online). Illustration of the edge locations of the v = 1/2 Moore-Read state and the v — 2 integer quantum 
Hall state. Although the positive background charge and total electron charge are 5 times that of the electrons forming the 
Moore-Read state, the edge is hiding behind the integer quantum Hall edge in the lowest Landau level. The Moore-Read edge 
may thus be protected since it is farther away from the electrostatic edge. 

where cj„ is the electron creation operator for the first excited Landau level (ILL) single electron state with angular 
momentum m. V^ ln 's are the corresponding matrix elements of Coulomb interaction for the symmetric gauge, and 
UmS the matrix elements of the confining potential. Additional details of this model can be found in Ref. [29| where 
we studied edge reconstruction of Abelian fractional quantum Hall states at different v, including explicit expressions 
of U and V and an illustration of the electrostatic configuration associated with He (see Fig. 1 of Ref. l29t ). 

The confining potential we use here is motivated by the <5-doping technique in 2DEG fabrication. For GaAs/AlGaAs 
heterostructures, silicon impurities are deposited in an atomically thick layer at a distance d ~ 1000 A above the 
interface where 2DEG is located to reduce impurity scattering. Therefore, we model the background potential arising 
from these ionized silicon impurities, which ensure charge neutrality in the samples. Even at an electrostatic level 
it is clear that d parameterizes the strength of the confining potential. At small d the potential is strong and also 
sharp near the edge, while at large d the potential is weak and smooth near the edge. Alternatively, one may tune the 
background charge density (right at the 2DEG plane) by smearing out the edge charge density as in an earlier study 
of edge reconstruction in integer quantum Hall liquids. 30 In the study of Abelian fractional quantum Hall liquids, 
we find that, e.g., the Laughlin-like state is stable up to d ~ 1.51b, beyond which edge reconstruction takes placed 
While we expect the parameter d appropriately characterizes the confining potential, we note the detailed sample 
structures and fabrication processes have an effect on how realistic the model is. 

To study the physics at v = 5/2, we explicitly keep the electronic states in the half-filled ILL only, while neglecting 
the spin up and down electrons in the lowest Landau level (OLL), assuming they are inert. The amount of positive 
background charge is chosen to be equal to that of the half-filled ILL, so the system is neutral. The choice of a disc 
radius R — y/iNl-Q, where N is the number of electrons in the ILL, guarantees that the disk encloses exactly 2N 
magnetic flux quanta, corresponding to v = 1/2 in the ILL. This is a simplification of the real system. In reality, the 
background charge equals the total electron charge of both the half-filled ILL and the filled OLL electrons. The latter 
neutralizes 4/5 of the background charge in the bulk, but this neutralization effect is incomplete near the edge due 
to finite d. Furthermore, the location of the OLL edge is different from that of the ILL electrons (see Fig. [IJ. The 
physical consequences of these effects will be discussed in Section IVIIIl In this study, we do not consider the finite 
thickness of the electron layer, which softens the Coulomb interaction and can be studied using the same numerical 
method, albeit time-consumingly. 

III. COMPETING GROUND STATES 

Taking advantage of the rotational invariance of the system, we diagonalize the Hamiltonian [Eq. |T])] for each 
Hilbert subspace with a total angular momentum M, and correspondingly obtain the ground state energy E(M). 
The global ground state is defined as the ground state with the lowest energy E(M gs ), whose corresponding angular 



5 



1.5 |T 




u — 1 — ^ — 1 ^ — 1 ^ — 1 — ^ 

0.0001 0.001 0.01 0.1 1 
X 



FIG. 2: (color online). Total angular momentum of the global ground state as a function of the mixing parameter A and 
background charge distance d for 12 electrons in 22 orbitals with the mixed Hamiltonian [Eq. ([T}]. The ground state in region 
(ii) has the same M ga = 126 as the Moore-Read (or Pfaffian) wave function. The ground state in regions (i) and (iii) are 
believed to be stripe phases. In region (iv), the ground state is a candidate for the so-called anti-Pfafhan state (see Sec. IVIll 
for detail). 

momentum is M gs . In our approach the ground state angular momentum is a result that comes out of the calculation, 
rather than a parameter fixed a priori based on the property of the state that one is interested in. Therefore, we can 
quantitatively analyze the stability of the ground state. 

Figure[2]is a phase diagram that shows the total angular momentum of the global ground state M gs for 12 electrons 
in 22 orbitals with the Hamiltonian in Eq. (fT]). We vary the mixing parameter A and the background charge distance 
d. The Moore-Read state for 12 electrons has Mmr = N(2N — 3)/2 = 126. In the small A limit, the ground state 
around d = 0.6-0.7 persist to have M gs — 126. To be precise, the ground state is stable for 0.51 < d < 0.76 for the pure 
Coulomb case A = 0j2& On the other hand, the range extends as A increases, since three-body interaction favors the 
Moore-Read state. The two regions with M gs = 121 and 136 surrounding the Moore- Read ground state are believed to 
be stripe phases. They can be represented by two strings of and l's \M gs = 121) = |1000001111111111100000) and 
\M gs = 136) = 1 1100000001 1111 11 11 1000), respectively. The and l's are the occupation numbers of single-electron 
angular momentum eigenstates (smaller angular momentum orbitals to the left). Alternatively, one can understand 
such a string as the Slater determinant of the corresponding single-electron angular momentum eigenstates labeled 
by 1. At this system size, numerical ground states have an overlap of about 30-40% with the corresponding Slater- 
determinant states in their range of stability. For very small d (d « 0.11b), M gs can jump to 110 for A < 0.01, which 
is believed to be a finite-size artifact. On the other hand, there is a region with ground state M gs = 146, which the 
authors already showed in Fig. lb of Ref. [26|. Toward the pure Coulomb case, this ground state is stable over a range 
of d twice as large as that for the Moore-Read state. We speculate this is related to the so-called anti-Pfaffian state 
discussed very recentlyj 17 ' 18 We will discuss this state in greater detail in Sec. IVIII 

In the absence of three-body interactions, the overlap between the ground state wave function and the Moore-Read 
wave function, \ (^ gs (M gs = 126)|\E'm_r) | 2 , is about 0.5 for the Coulomb interaction (it jumps up to 0.7 when we tunc 
the V\ pseudopotential 26 ). While this is quite substantial considering the already quite large size of Hilbert subspace, 
it is significantly smaller compared with the Laughlin state at v = 1/3 at comparable system size. Combined with the 
narrow window of d within which M gs — 126 in the pure Coulomb interaction case, these suggest that the Moore-Read 
state may be quite fragile when system parameters are varied. This is consistent with earlier numerical work on the 
torusi and the experimental observation that the FQH state at v — 5/2 disappears in a tilted magnetic field for 
modest tilting angle, even though the state is believed to be spin polarized. We note the phase boundaries in the 
small A limit persist in the small negative-A regime (—0.02 < A < 0). This suggests the Moore- Read-like ground state 
with pure Coulomb interaction is stable against a small attractive three-body interaction, which may arise, e.g., due 
to Landau level mixing. 
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IV. NON-ABELIAN QUASIHOLES WITH ELECTRIC CHARGE +e/4 IN THE MOORE-READ STATE 

Considering the relatively small overlap and rather narrow window of stability in d, one might wonder if the 
M gs = 126 ground state is indeed in the same universality class as the Moore-Read state. To answer this question 
we must study whether the elementary excitations of this state have the same properties as those of the Moore- Read 
state. In this section we study the the quasihole excitations of this state by introducing a local potential, possibly 
induced in experiments by the tip of an atomic force microscope, for example. The next section will be devoted to 
study of the edge excitations. 

As the ground state of the microscopic model is very sensitive to the parameters of the system, such as the 
background confining potential (by tuning d) and the weight (A) of the three-body interaction H^g 1 one may ask if 
additional features besides the total angular momentum can offer further support that the ground state is in the same 
phase as the Moore- Read state. In fact, one of the most striking properties of the Moore-Read state is that it supports 
charge ±e/4 quasihole/particle excitations, which carry half the charge of a Laughlin quasihole/particle at this filling 
factor. They obey non-Abelian statistics, and their existence implies that electrons are paired in the ground state (in 
the same way that observation of h/2e vortices indicate electrons are paired in superconductors). We note that the 
Halperin 331 stated also supports e/4 charge. But it is a bilayer state with 1/4 filling in each layer, thus e/4 charge 
is not as surprising, as one can get it by threading a flux quantum through one layer only. 

In an earlier study^ we have demonstrated that a short-range impurity potential at the origin Hy*/ = WcqCq can 
induce such a fractionally-charged quasihole, in the presence of some three-body potential. In a system of 12 electrons 
in 24 orbitals (as well as a smaller system of 10 electrons in 20 orbitals), we found for large enough W, a quasihole of 
charge +e/4 can appear at the origin. This is reflected in the depletion of 1/4 of an electron in the total occupation 
number of electrons at orbitals with small angular momenta, and in the change of ground state angular momentum 
from M gs = N(2N - 3) /2 to N(2N - 3)/2 + N/2, in agreement with that of the Moore-Read state with the quasihole 
located at the origin: 

*Mfl {Zl, Z 2 , Z N ) 

If W is increased further, a +e/2 quasihole (which is a Laughlin quasihole, equivalent to two +e/4 quasiholes^ 3 .) appears 
at the origin in the global ground state, whose total angular momentum further increases to N(2N — 3)/2 + N, in 
agreement with the variational wave function 

^~mr{zi,Z2,-,z n ) = ^11^ Vmr{zi,Z2,...,z n )- (13) 

In Fig. [3l we show the electron densities of the +e/4 quasihole and the corresponding Moore-Read ground state for 
30 electrons obtained by Monte-Carlo simulations. We note that the counterpart of Eq. [12] on the sphere 34 would 
represent two +e/4 quasiholes on the opposite poles of the sphere. In the following, we proceed to explore the existence 
of the +e/4 quasihole in a larger parameter space, including cases without three-body potential (A = 0). Since a +e/4 
quasihole with Abelian statistics can arise from a strong-pairing state (instead of the weak-pairing Moore- Read state), 
we will discuss the statistics of the quasihole in Sec. I VII after we discuss the edge excitations of the ground state. 
We attempt to trap a quasihole at the origin by introducing a Gaussian impurity potential: 

H w ,a = Wj2^p(-m 2 /2a 2 )clc m , (14) 

tn 

where a (multiplied by lg) is the range of the potential. Note H\y — WcqCo is the short-range limit (a — > 0) of 
the potential in Eq. 1141 Therefore, the additional parameter a allows us a more complete search. Two of us and 
a co-worker have also been studying the effect of the range and shape of the potential on the excitation of ±e/3 
quasiholes/quasiparticles in a Laughlin v — 1/3 liquid. 32 For a ~ 2.0, the weakest strength of the Gaussian potential 
that supports the quasihole state as the global ground state is found insensitive to the confining potential (or d in our 
model) <22, In the Moore-Read case, studies also suggest a ~ 2.0 is optimal for the generation of quasiholes, as in its 
vicinity the quasihole states can remain to be the global ground state even in the pure Coulomb case. 

Figure 2] shows the global ground state angular momentum as a function of the mixing parameter A and the tip 
potential strength W for 12 electrons in 22 orbitals, with the mixed Hamiltonian in Eq. |T]) and the Gaussian tip 
potential in Eq. (fT4")) . Here, we choose d = 0.71b and a = 2.0. To be specific, we expect the Moore- Read state with 0, 
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FIG. 3: (color online). The densities of the +e/4 Moore-Read quasihole state (Eg . I12|) and the Moore-Read condensate (Eq.|3]) 
for a 30-electron system. 

1, and 2 (+e/4) quasiholcs to have total angular momenta of 126, 132, and 138 respectively. For small W (W < 0.03), 
we find M gs = 126 for the global ground state, which is the same as the Moore-Read state. When there is enough 
three-body interaction (A > 0.025), as we increase W, M gs first jumps to 132 (W > 0.05), and then to 138 (W > 0.2) 
as W increases; this is exactly what one expects when the system first traps a single +e/4 quasihole, and then two 
+e/4 quasiholes. However, for smaller A, there is an additional region with M gs = 126 around W = 0.5, separating 
the one-quasihole and two-quasihole regions. This region turns out to be a stripe state, characterized by the binary 
string 10000011111111111100000). Careful analysis suggests that near A = 0.025 and W = 0.5, the energies of the 
three states with different total angular momenta are very close to each other and therefore extremely sensitive to the 
parameters. Despite this complication, we point out that the trapping of a single +e/4 quasihole by a local potential 
is a robust feature of the ground state, which persists to the pure Coulomb case along the lower boundary (W « 0.03), 
at least for finite potential width of a = 2.0. This strongly suggests that the ground state with M gs = 126 is indeed 
in the universality class of the Moore-Read state. 

We note that a 5-function trapping potential (<r — > in our Eq. (fT4")0 was used to generate quasiholes on a sphere 
by Toke et al£^ They were unable to isolate individual +e/4 quasiholes for either the pure Coulomb or pure three- 
body interactions, while in our earlier work2& we succeeded in doing that on a disc for some mixture of Coulomb and 
three-body interactions, using the same trapping potential. One advantage of disc geometry is that one can create a 
single +e/4 quasihole in the system, while on a sphere (or torus) such quasiholes must be created in pairs, and their 
interaction complicates the matter. Here we demonstrate that a single +e/4 quasihole can also be generated and 
isolated for pure Coulomb interaction, with some finite-range trapping potential. We have not, however, been able 
to do that with the 5-function trapping potential. This suggests that such quasiholes have relatively large size, and 
its trapping and manipulation will be sensitive to the details of the trapping potential. Thus experimentally one may 
need to optimize the trapping method in order to generate and manipulate them. 

V. EDGE EXCITATIONS OF THE MOORE-READ STATE AND THE INTERFERENCE 

EXPERIMENTS 



In addition to quasihole/quasiparticle properties, another way to probe the topological order of FQH liquids is to 
study their edge excitations, which are also of vital experimental importance. For comparison, the Laughlin state 
supports a single branch of bosonic chiral edge mode, whose properties have been studied in tunneling experiments. 4 
For the Moore-Read state, in addition to a bosonic mode whose properties are very similar to that of the Laughlin 
state, a neutral fermionic branch of excitations has been predicted^^ this fermionic branch is closely related to 
the non-Abelian nature of the state. The existence of both branches makes the low-energy excitation spectrum of a 
microscopic model at v = 5/2 richer and their experimental consequences more interesting^ 

In our earlier study^ we have observed both branches of modes for a mixed Hamiltonian, and demonstrated that 
a single +e/4 quasihole in the bulk changes the boundary condition of the fermionic mode, clearly indicating the 
non-Abelian nature of the quasihole. In this section we provide further details of the analysis of the spectra, and 
study how the spectra evolve as the interaction is varied, especially toward the pure Coulomb interaction. 
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FIG. 4: (color online). Ground state angular momentum as a function of trapping potential and three-body interaction strength. 
We have 12 electrons in 22 orbitals, with the mixed Hamiltonian [Eq. {l}] and the Gaussian tip potential [Eq. (|14|l ]. A and W 
characterize the three-body interaction and tip potential strength respectively. The background charge is fixed at d = 0.71b 
above the electron layer (the ground state has M gs = 126, same as the Moore-Read state in the absence of the Gaussian 
potential). For large enough A, as the tip potential strength W increases, states with M gs = 132 or M gs = 138, believed 
to contain a +e/4 quasihole or a +e/2 quasihole, become the global ground state. For small A, another ground state with 
M gs = 126 (which is a stripe state with occupation pattern |0000011111111111100000}), separates these two quasihole states. 



A. Edge spectrum of a Hamiltonian with mixed electron-electron interaction 

In this subsection, we demonstrate a clear separation of the fermionic and bosonic modes for the Moore- Read state, 
and try to obtain their velocities for A = 0.5. We will then try to extend the results to the pure Coulomb case in the 
next subsection. We begin by recalling the procedure to extract edge mode dispersion in the simpler Laughlin case 
at v = 1/3, where there is only one bosonic branch of edge mode. Then, we apply a similar analysis to the Pfaffian 
case, where we have a fermionic branch of edge mode in addition to a bosonic one. Of course, unlike the Laughlin 
case, here we need to rely on several reasonable assumptions, which can be justified post priori. 

In an earlier work^ we studied the energy spectrum of the electron system at v — 1/3, trying to identify the single 
bosonic branch predicted by the chiral Luttinger liquid theory. 1 The basic idea is that the low-lying excitations of the 
quantum Hall system at v = 1/3 can be described by a branch of single-boson edge states with angular momentum 
I (I = 1, 2, 3, ...) and energy e&(Z). Therefore, we can label each low-energy state by a set of (bosonic) occupation 
numbers {n(Z)}, whose total angular momentum is 

M = M + AM = M + n ( l ) l i ( 15 ) 

I 

and energy 

E = E + AE = E + J2 n ( l Ml), (16) 
i 

respectively, where Mo and Eq are total angular momentum and energy of the corresponding ground state. In 
Eq. (|16p we assumed the interactions between the excitations are negligible, which turns out to be an excellent 
approximation. Being edge excitations, such states can be independently verified by calculating the squared matrix 
elements T[{n(Z)}] = \(ip^ n my(N + ^)\cl N+AM \^o(N))\ 2 numerically in the microscopic model, and comparing them 
with the predictions of the chiral Luttinger liquid thcoryj 29 i 37 Note that Mq(N + 1) — Mq(N) = 3N is the difference in 
total angular momenta between the N- and (N + l)-electron ground states. As shown in Ref. [29l even in the presence 
of background confining potential, the ansatz of Eqs. (TT5)) and (fTo|) can be used to unambiguously identify the bosonic 
mode energies ej,(i), given that edge excitations are not significantly mixed with bulk excitations. The calculation of 
T[{n(Z)}], while not necessary, does ensure us the correct identification of these excitations as edge states. 

Encouraged by the success of identifying the edge mode dispersion and even predicting the energies of edge excita- 
tions in the Laughlin case, we apply the same analysis to the Moore-Read state. The complication is that, in addition 
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FIG. 5: (color online). Low-energy excitations AE(AM) from exact diagonalization for TV = 12 electrons in 26 orbitals in the 
ILL (corresponding to v = 1/2) for the mixed Hamiltonian in Eq. ((TJ with A = 0.5. The neutralizing background charge for 
the Coulomb part is deposited at d — 0.61b above the electron plane. The red solid bars, the black dashed bars, and the blue 
dotted bars mark fermionic, bosonic, and mixed edge excitations, respectively (see Appendix [X] for detail). Bulk excitations 
are represented by thin dotted bars (green). 



to the bosonic mode, we also have a fermionic mode, and thus the convolution of fermionic and bosonic excitations. 
Figure [5] shows the low-energy excitations AE(AM) for N — 12 electrons in 26 orbitals in the ILL for the mixed 
Hamiltonian with A = 0.5 and d — 0.61b- A gap at around AE = 0.1 is clearly separating the energy spectrum into 
a low-energy section and a higher-energy one. The numbers of the low-energy states for AM = M — Mq = 0, 1, 2, 
3, and 4 are 1, 1, 3, 5, and 10, respectively, agreeing perfectly with the numbers expected for the Moore-Read state 
by earlier theoretical workjiiii Notably, the lowest two levels for AM = 4 lie very close to each other. Based on the 
agreement in numbers, we are tempted to call them edge states; but further confirmation comes from the separation 
and identification of bosonic mode and fermionic mode, as we show below. 

We assume each low-energy excitation can be labeled by two sets of occupation numbers {nb(lb)} and {n>f(lf)} for 
bosonic and fermionic modes with angular momenta lb, If, and energies e&, e/, respectively. rib(l b ) are non- negative 
integers while rif(lf) = 0,1. Since the fermionic edge excitations are Majorana fermions that obey antiperiodic 
boundary conditions^ If must be positive half integers, while for bosonic mode l b are integers. In addition, the 
total fermion occupation number n f(lf) f° r eacn state must be an even integer because each fermionic excitation 
contains an even number of Majorana fermion modes due to their pairing nature. The angular momentum and energy 
of the state, measured relatively from those of the ground state, are 

AM = J2Mh)h + Y, n f(h)lf, (17) 

l b If 

AE = ^2n b (l b )e b (l b ) + ^2nf(lf)€f{l f ). (18) 

h h 

The details of the analysis on the data of Fig. (O are presented in Appendix [A] Here we summarize the results 
in Tableland Fig. [Ha). Interestingly, the fermionic dispersion curve is monotonic and can be well fit by a straight 
line passing the origin, allowing us to obtain the neutral fermionic velocity v n = dtf/dk ~ 0.0016i?e 2 /(e/s?i), where 
the disc radius R — 2^/NIb, and we have the conversion from the angular momentum to the linear momentum along 
the edge fc = AM/R. For typical GaAs systems, we obtain »„Ri2x 10 5 cm/s. On the other hand, in contrast to the 
roughly linear dispersion of the fermionic branch, the energy of the bosonic branch bends down (despite a much bigger 
initial slope or higher velocity), suggesting a potential vulnerability to edge reconstruction in the bosonic branc h 29 ' 38 . 
This is not surprising since the bosonic mode is charged; as a result its velocity is dominated by the long-range 
nature of the Coulomb interaction in the long- wavelength limit, but at the same time it is also more sensitive to the 
competition between Coulomb interaction and confining potential which can lead to instability at shorter wavelengths. 
If we assume the curve is linear for k < 1/R, we can estimate v c « 3 x 10 6 cm/s for GaAs. 

Using these e&'s and e/'s (a total of 8 energies), we can re-construct the whole low-energy spectrum of the system 
up to AM = 4 (a total of 20 states), which agrees well with the actual spectrum 2 ^ (in fact, we can extend the 
construction to AM = 5 and obtain very satisfactory agreement for most states, which do not involve edge modes 
with larger momentum). The consistency justifies our analysis based on the assumption that the interactions between 
excitations are negligible and further supports our central result in this section, namely the fermionic mode is well 
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tb(h) 
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£/('/) 


1 


0.022659 


1/2 


0.000324 


2 


0.030057 


3/2 


0.001676 


3 


0.029908 


5/2 


0.004117 


4 


0.024668 


7/2 


0.006011 



TABLE I: Dispersion energies of both bosonic and fermionic modes at small momenta for TV = 12 electrons at half filling 
(in 26 orbitals) in the first Landau level. The system has a Hamiltonian of 50% Coulomb interaction and 50% three-body 
interaction. The background charge is placed at d — 0.6. Based on these energies, we can construct the complete low-energy 
(edge) spectrum of the 12-electron Pfaffian state up to AM = 4 (Ref. 26). 




1 2 3 4 0.5 

AM X 



FIG. 6: (color online), (a) Dispersion curves of bosonic (et) and fermionic modes (e/) for N = 12 electrons at half filling (in 
26 orbitals) in the first Landau level for A = 0.5. These energies can be used to construct the complete edge spectrum for the 
12-electron system up to AM = 4.— (b) Dispersion curves for the same system with A = 0.1, or less three-body interaction, 
(c) Bosonic (v c ) and fermionic velocities (v n ) extrapolated to the pure Coulomb case. We obtain the A = 0.5 and 0.1 points 
from fitting the fermionic modes to a straight line in (a) and (b), respectively. At A = 1 (pure three body interaction, v n = as 
all edge states have zero energies. We thus obtain v c = 0.046 and v n = 0.0036, in units of (i?e 2 )/(e/_g?i), for the pure Coulomb 
case. 

separated from the bosonic mode and has much lower energy and velocity. 

B. Toward the pure Coulomb interaction 

The ultimate goal of our work is, of course, to understand the low-energy spectrum with pure Coulomb interaction, 
or at least with less three-body interaction. Looking at the energy spectra for A = 0.1 (Fig. [7]) and 0.0 (Fig. [8]), we 
fail to observe a gap separating edge and bulk states, as in Fig. [5] One interesting question is, as they start to have 
similar energies, whether bulk states and edge states are mixed. However, without the gap, it is difficult to identify 
each eigenstate as a bulk state, a specific edge state, or a mixture of edge and bulk states. To allow the identification, 
we calculate the overlaps between the eigenstates for A = 0.5, which we have already analyzed, and the eigenstates 
for A = 0.1 and 0.0. We assume the eigenstates evolve smoothly as A decreases, which turns out to be the case as our 
analysis will show. Thus, we can trace the edge states identified for A = 0.5 and sort them out from all eigenstates in 
the pure Coulomb case by calculating overlaps; the sorting is otherwise impossible. In particular, we are interested 
in the evolution of fermionic edge states, which play an important role in understanding the non-Abelian nature of 
the Moore-Read state. 

We leave the details of the approach to Appendix IBl but highlight the main results here. We first look at A = 0.1. 
Figure [7] shows the low-energy excitations for 12 electrons in 26 orbitals in the ILL for the mixed Hamiltonian in 
Eq. with A = 0.1. The neutralizing background charge for the Coulomb part is deposited at d = 0.6/b above the 
electron plane, just as in the previous case. We find the fermionic edge excitations (red solid bars) are well separated 
from bulk and other edge excitations, as there is clearly a spectral gap around AE = 0.02. Another observation is for 
AM = 4. Here the two fermionic excitations can be significantly mixed with each other, as the two states for A = 0.1 
have roughly equal overlap with the two for A = 0.5. 
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FIG. 7: (color online). Low-energy excitations AE(AM) from exact diagonalization (solid lines) for TV = 12 electrons in 26 
orbitals in the ILL (corresponding to v = 1/2) for the mixed Hamiltonian in Eq. |T} with A = 0.1. The neutralizing background 
charge for the Coulomb part is deposited at d = 0.6Zb above the electron plane. The red solid bars, the black dashed bars, 
and the blue dotted bars mark fermionic, bosonic, and mixed edge excitations, respectively. Bulk excitations are represented 
by thin dotted bars (green). While the bosonic edge excitations mix significantly with the bulk excitations, the fermionic edge 
excitations are still well separated from the rest (see Appendix IB II for detail). 
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FIG. 8: (color online). Low-energy excitations AE(AM) from exact diagonalization (solid lines) for N = 12 electrons in 26 
orbitals in the ILL (corresponding to v — 1/2) for pure Coulomb interaction (A = 0.0). The neutralizing background charge 
for the Coulomb part is deposited at d = 0.6/b above the electron plane. The red solid bars, the black dashed bars, and the 
blue dotted bars mark fermionic, bosonic, and mixed edge excitations, respectively. Bulk excitations are represented by thin 
dotted bars (green). In this case, fermionic edge excitations also mix with the bulk excitations (see Appendix lB 2l for detail). 



Similar to what we have done in the previous subsection, we can extract the bosonic and fermionic mode energies 
for A = 0.1, plotted in Fig. [6fb). Interestingly, the figure looks exactly like that for A = 0.5, except the energy scales 
are roughly doubled for a higher percentage of Coulomb interaction. The bosonic curve bends down slightly further. 
In this case, v n (\ = 0) = 0.0033i?e 2 /(^s?i)- 

Figure [5] shows the low-energy excitations for pure Coulomb interaction and the confining potential with d = 0.6 
for 12 electrons in 26 orbitals. Unfortunately, there is no clear distinction between edge modes and bulk modes at this 
system size. The situation here is similar to a related study on a rotating Bose gas^ After calculating the overlaps, 
we find the energies of the lowest fermionic edge excitations are around 0.02, but there are lower eigenstates which 
originate from bulk excitations. We note that a recent DMRG study suggest the bulk gap of the fractional quantum 
Hall liquid at v = 5/2 is approximately 0.034°- 

Unfortunately, we can no longer extract meaningful results for the bosonic and fermionic dispersion curves for the 
pure Coulomb interaction as we did for A = 0.5 or 0.1. We believe the reason is that the fermionic mode starts to mix 
with bulk states, distorting the dispersion curves (see Appendix IBl for detail). Nevertheless, we can extrapolate the 
bosonic and fermionic velocities from the two finite-A values, along with the fact that the velocities are zero for the 
pure three-body case A = 1, because all edge states have zero energy in iJ 3B 4i This also suggests that the velocities 
should be roughly proportional to the weight of Coulomb interaction 1 — A, which is indeed what we find in Fig.^c). 
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The extrapolations give v c = 0.046 and v n — 0.0036 for the pure Coulomb case (A = 0), in units of Re 2 /(elsK). In 
GaAs systems they are v c « 5 x 10 6 cm/s and v n w 4 x 10 5 cm/s, respectively. We can check the validity of the 
numbers using Fig. [8j for example, the bosonic state at AM = 1 is at E = 0.0476, very close to 0.046 based on the 
value of v c . Similar comparisons also find v n w .004(i?e 2 ) / '(el bK) is in reasonable agreement with the energies at 
AM = 2 and 3. 

We close this subsection by noting that the mixing between bulk and edge excitations seen here for pure Coulomb 
interaction or small A is a finite size effect. The edge excitations, which are gapless in the thermodynamic limit, 
have a finite gap due to the existence of a minimum momentum k dictated by system size. When this gap is larger 
than the bulk excitation gap (which is quite small for the Moore- Read-like state with pure Coulomb interaction), 
mixing between the two types of excitations occurs. They will ultimately separate as system size increases, in the 
long-wave length limit k — > 0. Adding the three-body interaction makes this separation occur at smaller system size, 
by increasing the bulk gap without affecting the edge excitation energy much. Thus the effect of adding three-body 
interaction is similar to increasing system size, which allows us to extract useful information within the accessible 
system sizes. 



C. Implications on interference experiments 



Our numerical calculation suggests that the neutral mode velocity v n is much smaller than the charge velocity 
v c . A similar conclusion has been reached in an effective edge theory study, which also suggests the neutral velocity 
has a dynamic origin 4i The situation is somewhat similar to what happens in a ID Luttinger liquid of spin- 1/2 
electrons, where the velocity of the spin mode is in general lower than that of the charge mode, leading to the so 
called spin-charge separation. Here we coin a similar term, "Bose- Fermi separation" , to describe the separation of the 
velocities of charged bosonic and neutral fermionic edge excitations of the Moore-Read edge. 

In a Luttinger liquid, spin-charge separation is a main source of the decoherence of a single electron 42 . Physically 
this is because an electron carries both spin and charge; once it enters the Luttinger liquid however, its spin and 
charge components propagate with different velocities, leading to physical separation between the two after some 
decoherence time, and loss of integrity of the electron. 

The same physics is relevant to the fate of a charge ±e/4 quasihole/quasiparticle when it is propagating along 
the Moore-Read edge. A charge ±e/4 quasihole/quasiparticle carries both a bosonic component and a ferminic 
component; the former carries its charge while the latter is responsible for its non-Abelian nature. Similar to the 
case of an electron in a Luttinger liquid, we expect Bose-Fermi separation to be a main source of decoherence of such 
a non-Abelian quasihole/quasiparticle when it propagates at the edge. This raises a concern that such decoherence 
may destroy the interference pattern coming from the interference between charge ±e/4 quasiholes/quasiparticles in 
intcrferometry experiments proposed recentl y 21 ! 22 In a very recent work^ it was found that the decoherence length 
is indeed very sensitive to the velocities: 

1 fl/8 1/8V 1 
L <t> = TT^ — + — • 19 

As a result, in a double point-contact interferometer, the oscillatory tunneling current due to interference of ±e/4 
quasiholes/quasiparticles decays like I cx e -i / L *, where L is the distance between the two point contacts. It is clear 
from the equation above that is controlled by v n , when v n <C v c , and smaller v n leads to shorter L^. 

Based on our numerical results, we can estimate the constraints on the interferometry experiments due to deco- 
herence. In the pure Coulomb case, we use the bosonic and fermionic velocities extrapolated in Fig. [He). Assuming 
the experiments are done at a temperature of 10 mK and a magnetic field of 5 T,— we estimate ss 4 /im; this 
raises concerns on the appropriate inter-point-contact distance L in interference experiments. In fact, this may be a 
(perhaps overly) optimistic estimate, as we have not considered the errors due to finite system size and other realistic 
issues like filled lowest Landau level. Most importantly, the confining potential we use in our model (with parameter 
d/ls ~ 1) is much stronger than that for real systems ; 29 ' 38 real samples have much bigger d/ls, resulting in weaker 
confinement and thus smaller v n , leading to a smaller (see next section for further discussion on this point). Thus 
our estimate using parameters extracted from the specific model we use is best viewed as an upper bound of L^. 
Further investigation on this is thus needed. 

We close this subsection by noting that while Bose-Fermi separation has important consequences on the decoherence 
of charge ±e/4 quasiholes/quasiparticles, it does not affect charge ±e/2 Laughlin quasiholes/quasiparticles that 
only carry the Bose component. The interference pattern due to these Laughlin ±e/2 quasiholes/quasiparticles, 
unfortunately, does not exhibit the exciting non-Abelian behavior. Thus in interference experiments it is possible that 
while the interference due to charge ±e/4 quasiholes/quasiparticles is lost due to decoherence, one can still observe 
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FIG. 9: (color online). Low-energy spectra in a system of 12 electrons in 24 orbitals for the mixed Hamiltonian with A = 0.1 
and d — 0.71b- (a) In the absence of the external tip potential, the ground state (M = 126) is Moore- Read-like. There are 
0, 1, 1, 2 low-energy (fermionic) excitations for AM = 1-4. (b) In the presence of a Gaussian tip potential (W = 0.1 and 
a — 2.0), the +e/4 quasihole state emerge at M = 132. There are 1, 1, 2, 2 low-energy (fermionic) excitations for AA/ = 1-4. 
This suggests that a single +e/4 quasihole changes the fermionic mode spectrum. (The states of interest are marked by red 
solid bars.) 



an interference pattern due to charge ±e/2 quasiholes/quasiparticles, which is similar to that in Laughlin states. 
In addition to Laughlin quasiholes/quasiparticles, there are also charge ±e/2 quasiholes/quasiparticles that carry a 
neutral fermion (ip) but are also Abelian. Bose-Fermi separation does affect their propagation and thus suppresses their 
interference. Also the added fermion component makes tunneling of such ±e/2 quasiholes/quasiparticles irrelevant^ 
further reducing their importance. 



VI. NON- ABELIAN NATURE OF +e/4 QUASIHOLE AND POSSIBLE INSTABILITY OF FERMIONIC 

MODE AT THE EDGE 



Now armed with the capability of exciting quasiholcs as well as the knowledge of edge modes, we are in a position 
to reveal the non- Abelian nature of a +e/4 quasihole by studying the change of fermionic edge states in the presence 
of the quasihole. Such a change has been reported in an earlier paper by three of the authors 26 for A = 0.5, d = 0.51b 
and a short-range tip potential. Here we are presenting a case with less three-body interaction (A = 0.1), weaker 
confinement (d = 0.71b), and a Gaussian tip potential. 

Figure [9] shows the low-energy spectra in a system of 12 electrons in 24 orbitals for the mixed Hamiltonian in Eq. ((T|) 
with A = 0.1 and d = 0.71b- In the absence of the external tip potential (Fig. [5^), the ground state (M — 126) is 
Moore- Read-like, as we can also read from the phase diagram in Fig. [5] (albeit in 22 orbitals). The excitation spectrum 
clearly has a gap up to roughly 0.03, consistent with the result from numerical DMRG calculations^ Inside the gap, 
there are 0, 1, 1, 2 low-energy excitations for AM = 1-4 (marked by red solid bars); the numbers agree precisely with 
the number of fermionic states, as discussed in Sec. [V] Comparison with Fig. [7] suggests that the fermionic mode 
dispersion gets distorted by the increased d (smoother confinement). 

In the presence of a Gaussian tip potential (W = 0.1 and a = 2.0), a new ground state emerges at M = 132, 
reflecting the fact that a +e/4 quasihole has been trapped by the tip potential. Now there are 1, 1, 2, 2 low-energy 
(fermionic) excitations for AM = 1-4. For any M, there is an energy gap of at least 0.016 separating the fermionic 
edge states and the rest. The results suggest that a single +e/4 quasihole changes the fermionic mode spectrum, a 
remarkable feature due to the non- Abelian nature of the +e/4 quasihole. More precisely, a quasihole carries a a field 
of the Ising conformal field theory, and its presence changes the boundary condition of the edge Majorana fermion 
mode from being antiperiodic to periodic. This leads to a shift of the angular momentum quantum numbers of the 
lowest energy fermionic edge excitations Since the non- Abelian properties of such quasiholes are exclusively due to 
the a degree of freedom it carries, observing such a change of boundary condition directly confirms the non- Abelian 
nature of the quasihole. 

We note that for this particular set of parameters, two of the low-energy femionic edge states (AM = 1 and 3) 
actually have very small negative energies measured from the single quasihole state with M = 132. This is sensitive to 
the choice of parameters; for sharper confinement (with a cutoff of 22 orbitals), the quasihole state is the global ground 
state. This does suggest that there are potential instabilities in the fermionic mode; such possible instabilities and 
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their consequences remain to be investigated. We note that even if such instabilities do not occur, the closeness of the 
fermionic excited state energies to zero suggests the neutral velocity can be even smaller with smoother confinement, 
which can jeopardize the estimate we made in the previous section on the dephasing length relevant to double point- 
contact interference experiments. 

VII. DISCUSSION ON THE POSSIBLE ANTI-PFAFFIAN STATE 

In Sec. IIII1 we mentioned a stable ground state with the same angular momentum quantum number as the recently 
proposed anti-Pfaffian state.— ^ In this section, we discuss how this quantum number is determined. 

Suppose we have a system of N electrons in N or t orbitals. The Pfaffian state, or the Moore-Read state, has a total 
angular momentum of N(2N — 3)/2. In order to be able to accomodate this state, we need Norb > 2N — 2. From 
another angle, we can equivalently view the system as having Nh = N orb — N holes in N orb orbitals. The simplest 
version of the anti-Pfaffian state, by definition, is the Pfaffian state formed by all the holes; this is possible as long as 
N orb > 2iVft — 2, or N orb < 2N + 2. The total angular momentum (for holes) is thus 

M h = ^ (2N h - 3) = Norb ~ N [2{N orb — N) — 3]. (20) 

The total angular momentum in the original electron basis is 

M A p = N orb {N orb - l)/2 - M h , (21) 

where the first term is the contribution from the electron background that fully occupies all N oro orbitals, and the hole 
contribution Mh is negative because a hole removes an electron from an occupied orbital. For N = 12 and N orb = 22, 
we find Nh = 10, Mh = 85, and Map = 146. This is exactly the total angular momentum of the ground state in the 
region (iv) in Fig. O Furthermore we found that increasing the three-body interaction enhances the Pfaffian state 
and suppresses this state; this is consistent with our speculation that this is the anti-Pfaffian state. 

This is, however, not definitive evidence, as there are competing states with the same quantum number. In 
particular, a 12-electron stripe-like state represented by the binary string |1100000000111111111100) has the same 
angular momentum 146 and very low energy. Analysis of the system with N = 12 electrons in 22 orbitals with pure 
Coulomb interaction reveals a large overlap (0.35) between the numerical ground state and the stripe state. 

We now explore more general possibilities of the anti-Pfaffian state by increasing N orb from 22 to 24. If the same 
10-hole anti-Pfaffian state were to be realized, the two extra holes would be at the two outermost orbitals, and the 
ground state will have the same quantum number. Stripe or other states, on the other hand, are more likely to respond 
to the change of boundary. In our numerical calculation, we indeed find the global ground state still has M to t = 146 
for d — 1.2. In addition, the overlap between the ground state and the stripe state discussed above decreases to 0.13. 
This seems to suggest that the stripe phase is favored by sharp (hard- wall) confinement, as the two outermost orbitals 
are unoccupied. With smoother confinement, a different state, which we speculate is related to the anti-Pfaffian state, 
emerges. At d = 1.5, the ground state momentum increases to 151; this can not be easily explained by a simple stripe 
phase. On the other hand, it could be explained as the anti-Pfaffian state with one +e/4 quasihole. 

These findings suggest that, due to the presence of a confining potential, the Pfaffian and anti-Pfaffian states 
have different energies, and may be realized at different confining potential strengths, without invoking effects of 
Landau level mixing. This is because the confining potential breaks the particle-hole symmetry. Based on our model 
calculation, we speculate that the Pfaffian state is stable as the ground state for strong confinement (sharp edge and 
smaller d), while the anti-Pfaffian state is stable for weak confinement (smooth edge and larger d). It is worth pointing 
out that the two phases are separated and strongly influenced by a stripe phase. Whether this is a generic feature 
or a finite-size artifact cannot be resolved in the current work. Since the anti-Pfaffian is stable around d — 1.5/b, it 
opens another interesting possibility that edge reconstruction 38 may play a more important role in the anti-Pfaffian 
state. 



VIII. CONCLUDING REMARKS 

To summarize, we have studied a microscopic model of fractional quantum Hall liquids at filling fraction v = 5/2. 
The interaction between electrons are interpolated continuously between the limits of purely the three-body interaction 
and purely the Coulomb interaction. Another parameter we vary in our study is the strength of confinement potential, 
parameterized by the distance d separating the positive neutralizing background charge and the 2D electron gas layer. 
This enables us to reveal the nature of ground states and elementary excitations in the pure Coulomb interaction 
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limit and with semi-realistic confining potential. In particular, we find a Moore-Read-like state is realized in a small 
window of parameter space, with predicted properties. 

The Moore-Read-like ground state has an edge spectrum consistent with that of a charged bosonic mode and a 
neutral fermionic mode. The fcrmionic mode has much lower energy than the bosonic mode, implying the neutral 
velocity is at least an order of magnitude smaller than the charge velocity. This leads to a constraint on the dephasing 
length for charge ±e/4 quasiholes/quasiparticles: < 4 /im for typical experimental parameters, at T — 10 mK. 
This length is of crucial importance in double point-contact interference experiments. 

A local potential with a finite width (~ 21b, or about 20 nm), modeling an atomic force microscope tip, can induce 
exactly one charge +e/4 quasihole or one charge +e/2 (equivalent to 2 charge +e/4 quasiholes). From the change of 
fermionic edge mode when a single charge +e/4 quasihole is excited, we confirm the non-Abelian nature tof charge 
+e/4 quasihole. 

A ground state with the same quantum number as the recently proposed anti-Pfaffian state is stable in a weak 
and smooth edge confining potential. The state is found to be separated from the Moore-Read Pfaffian state by a 
stripe-like state in finite-size calculation. 

In the present work, we have used a semi-realistic model for the numerical calculations and attempted to obtain 
concrete numbers in experimental units, although further improvement is certainly possible and probably necessary. 
Among the effects we have neglected here, perhaps the most important is the presence of the electrons occupying the 
lowest Landau level (OLL), and their associated edge (see Fig. [T|). These OLL electrons have two effects ethat are not 
included in our study. The first is that the background charge needs to be equal to the total electron charge, not just 
those in ILL. While the additional charge is neutralized by the OLL electron charge in the bulk, this neutralization is 
incomplete at the edge, which results in a fringe electric fields that tends to destabilize the ILL edge through edge 
reconstruction. On the other hand, due to the cyclotron gap between the OLL and ILL, the ILL edge "hides behind" 
the OLL edge, and gets protected from instabilities by the OLL edge. Thus these two effects that we neglected impact 
the ILL edge in opposite ways, and further studies are needed to resolve which effect dominates, and the ultimate 
fate of the ILL edge. 

Nevertheless we do believe the numbers obtained from the present work can be of use as guidance to experimentalists 
who are interested in engineering samples and devices or in manipulating individual non-Abelian quasiholes in these 
devices. The parameters considered here, which describes the smoothness of the edge (related to the location of 
(5-doping in realistic epitaxially grown samples) and the size of an atomic force microscope tip are intimately relevant 
to experiments. For example, a momentum-resolved magneto-tunneling study has found that an epitaxially overgrown 
cleaved edge can realize the sharp edge limits With these realistic issues in mind, this work supports the possibility 
of topological quantum computing 45 - using fractional quantum Hall states, although the road ahead needs further 
exploration. 

An immediate follow-up study, which can strengthen the confirmations found in this work, is the study of the 
effects of the electron layer thickness, currently under exploration. A recent study by Peterson and Das Sarma 4e 
claims that finite layer thickness enhances the Moore-Read state, using the criterion of wave function overlap. It 
would be interesting to study the layer thickness effects in our more sophisticated model using criteria involving 
ground state energy, bulk and edge excitations. In addition, one also desires to look at the results in larger systems, 
where finite-size effects are weaker. Techniques to reduce the size of the Hilbert space using various truncation schemes 
are under development. 
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APPENDIX A: ANALYSIS OF EDGE EXCITATIONS AT A = 0.5 

In this appendix, we discuss in detail the analysis of the edge excitations for the mixed Hamiltonian with A = 0.5, 
plotted in Fig. [5] Below E = 0.1, we have 1, 1, 3, 5, and 10 states with respective angular momenta AM = M — Mq = 
0, 1, 2, 3, 4, which are well separated from the rest; we identify them as low-energy excitations below the bulk excitation 
gap. The sequence of numbers are those expected from edge excitations made of a chiral bosonic branch and a chiral 
fermionic branch. Therefore, we want to associate each of the 20 states with two sets of occupation numbers {rib(lb)} 
and {n-f(lf)} for bosonic and fermionic modes with angular momenta If,, If, and energies et,, ef, respectively. 
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TABLE II: Overlap matrix of the two systems with A = 0.5 (row) and 0.1 (column) for N = 12, N or b = 26, M = 127, and 
d — 0.6. The largest overlap between eigenstates for A = 0.1 and the lowest state (edge state) for A = 0.5 comes from the 
second lowest state for A = 0.1, with a value of 0.827 (underlined). 

Besides the ground state, it is not difficult to identify the only low-lying state at AM = 1 as the bosonic mode 
with energy £b(l) = AE(AM = 1) = 0.022659. We can thus identify all edge states at energies AE = n&(l)ei,(l) with 
corresponding momenta AM — n&(l). 

For AM = 2, we associate the highest-energy state with AE as 2et>(l)- There are two more states left, with energies 
ef>(2) and e/(l/2) + e/(3/2). There are thus two choices. However, given e&(l) ~ 0.02, it is reasonable to assume 
£f,(2) = 0.030057 is the higher one of the two. As a result, the fermionic state with the smallest momentum has much 
lower energy than the bosonic ones. Counting the energy states with nearly zero energy (or to be more precise, with 
AE < 0.01), we find 0, 1, 1, and 2 states for AM = 1, 2, 3, and 4, respectively. These numbers agree perfectly with 
the results expected for a single branch of Majorana fermion modeJ^- We thus assume these energies are sums of 
two Majorana fermion energies. For AM = 2, for example, we have already assumed AE = e/(l/2) + ei (3/2) for the 
only state. 

For AM = 3, we have five states. We continue to assume that the lowest one is a purely fermionic state with 
AE = ey(l/2) + e/(5/2). We can also identify two bosonic states with energies AE = 3ef,(l) ~ 0.066 and e&(l) + 
6f,(2) « 0.052. We also find one more from the convolution of both bosonic and fermionic modes with energy AE = 
e b (l) + 6/(1/2) + e/(3/2) « 0.024. The edge state left should then be the bosonic state with AE = e 6 (3) = 0.029908. 

The situation becomes more complicated at AM = 4, where we have two fermionic, five bosonic, and three 
convoluted edge excitations. It is easy to identify the convoluted excitations first, at energies e&(l) +e/(l/2) + e/(5/2), 
e&(2) + ef (1/2) + e/(3/2), and 2eb(l) + e/(l/2) + e/(3/2). Using ej,(i) for I = 1-3 obtained above, we can identify four 
bosonic states at energies 4e&(l), 2e&(l) + et>(2), 2e&(2), and e^l) + 6^(3). The only state with energy AE > 0.01 
is thus the remaining bosonic state with Cb (4) = 0.024668. Once again, the two fermionic states have much smaller 
energy AE — e/(l/2) + e/(7/2) and e/(3/2) + e/(5/2). We note that in order to write down a variational wave function 
for a pair of Majorana- Weyl fermions with momenta I > k, we need at least (2N + 1 — 1) orbitals. Therefore, by 
reducing the Hilbert space by using fewer orbitals, the hard- wall edge confinement will increase some fermionic mode 
energies, but leave others intact. This is a test that can unambiguously distinguish the two states. In particular, 
e/(7/2) will suffer from an energy increase when we reduce the total number of orbitals to 25, while e/(5/2) will 
remain roughly unchanged unless we further reduce the orbital number to 24 and below. We have observed this 
confinement effect in numerical calculations, which suggests the state with energy e/(3/2) + e/(5/2) is the lower of the 
two. This energy, together with the two fermionic excitations at smaller momenta, allow us to solve e/(l/2), e/(3/2), 
and e/(5/2). Consequently, the energy of the other state [e/(l/2) +ef(7/2)] allows us to solve for e/(7/2). The results 
are summarized in Table U and Fig.[5J 

Interestingly, the fermionic dispersion curve is smooth, monotonic and can be well fit by a straight line passing the 
origin, allowing us to obtain the neutral fermionic velocity v n . The bosonic dispersion curve, on the other hand, is 
non-monotonic and bends downward, which indicates a tendency toward edge reconstruction. 38 A very similar analysis 
can be perfromed for A = 0.1, which we leave out for brevity. The results are compared with A = 0.5 in Fig. [51 

APPENDIX B: NUMERICAL ANALYSIS OF THE EVOLUTION OF EDGE EXCITATIONS 

1. A = 0.1 

We discuss here in detail how we identify the edge excitations in Fig.[7]by calculating the overlaps between eigenstates 
for different A. Let us start with the simplest nontrivial case AM — 1 (Table [TO . The lowest excitation (state #1) 
for A = 0.5 has the largest overlap (0.827, underlined) with the second excitation (state #2) for A = 0.1. Meanwhile, 
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TABLE III: Overlap matrix of A = 0.5 (row) and 0.1 (column) for N = 12, N orb = 26, M = 128, and d = 0.6. The underlined 
numbers are the overlap between an edge state for A = 0.5 and the (likely) corresponding edge state for A = 0.1. 
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TABLE IV: Overlap matrix of A = 0.5 (row) and 0.1 (column) for N = 12, N orb = 26, M = 129, and d = 0.6. The underlined 
numbers are the overlap between an edge state for A = 0.5 and the (likely) corresponding edge state for A = 0.1. 

state #1 for A = 0.1 has large overlaps with states #2 and #3 for A = 0.5, both bulk excited states. Therefore, we 
can identify (AM, AS) = (1,0.0415) (state #2 for AM = 1) as an edge state for A = 0.1, with an overlap of 0.827 
with the corresponding edge state #1 for A = 0.5. 

For AM = 2, we find that states #1, #4, and #12 have significant overlaps with the lowest three edge states 
for for A = 0.5, as listed in Table [TTTJ As in Fig. A£(#12,AM = 2) = 0.0828 is roughly the sum of two 
AE(#2, AM = 1) = 0.0415. This simple addition law resembles the one we have found for edge states at one of the 
Laughlin filling fractions v = 1/3,— reflecting the conservation of energy and angular momentum. 

We may assume, based on this observation, that AE(AM — 3) = 3A£?(#2, AM = 1) is an edge state (not plotted 
in Fig. [7]). The other four edge states for AM = 3 are found, by comparing overlaps, to be states #1, #3, #4, and #20, 
according to Table|IV] We note the approximate equalities A£(#20, AM = 3) « AE(#2, AM = 1) + A£(#4, AM = 
2) and A£(#3, AM = 3) « AE(#2, AM = 1) + AE(#1, AM = 2). 

Similarly, we can identify six edge excited states in the lowest 20 eigenstates we have calculated for AM = 4. In 
addition, we can postulate the existence of another four edge states with excitation energies of 4A_E(#2, AM = 1), 
2A£(#2, AM = 1)+A£(#4, AM = 2), 2A£(#4, AM = 2), and AE{#2, AM = 1)+A£(#4, AM = 4), respectively. 
Again, the simple conservation law seems to work fairly well. We point out that the two fermionic edge states (#1 
and #2), whose energies are close, mix significantly with each other with respect to the A = 0.5 case. To a lesser 
extent, another two states (#3 and #4) also mix with each other. 
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TABLE V: Overlap matrix of A = 0.5 (row) and 0.1 (column) for N = 12, N orb = 26, M = 130, and d = 0.6. The underlined 
numbers are the overlap between an edge state for A = 0.5 and the (likely) corresponding edge state for A = 0.1. 
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TABLE VI: Overlap matrix of the two systems with A = 0.5 (row) and 0.0 (column) for N = 12, N orb = 26, M = 127, and 
d — 0.6. The underlined element of 0.403 is the overlap of the lowest state (edge state) for A = 0.5 (mixed system) and the 
seventh lowest state for A = 0.0 (pure Coulomb system). 
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TABLE VII: Overlap matrix of A = 0.5 (row) and 0.0 (column) for N = 12, N orb = 26, M = 128, and d = 0.6. The underlined 
numbers are the overlap between an edge state for A = 0.5 and the (likely) corresponding edge state for the pure Coulomb case. 

2. The pure Coulomb case 

We now move to the pure Coulomb case with A = and look for the eigenstates with significant overlap with 
the edge states in the A = 0.5 system. Surprisingly for AM = 1, the 7th lowest state has the largest overlap with 
the bosonic eigenstate in the corresponding subspacc for A = 0.5 (Table IVT1) . The overlap 0.403 is far from unity, 
but comparable to that between the Pfaffian state and the ground state of the pure Coulomb system. The six lower 
eigenstates (with the notable exception of the state #5), which have negligible overlaps with the edge state, have 
nonetheless significant overlaps with the lowest energy bulk excited states 2-5 for A = 0.5, indicating their bulk nature. 
The complexity of the Coulomb case is thus evident even for the AM — 1 case. 

The attempt to find all edge states, even with the overlap matrix calculation, is challenged by the following two 
difficulties. First, the edge states now have very large excitation energies, and thus a lot more eigenstates are needed 
for the search. Second, the overlaps with eigenstates for A = 0.5 fail to exhibit a clear one-to-one correspondence. 
In many cases, an edge states for A = 0.5 can have comparable overlaps with two eigenstates for A = 0, making the 
identification ambiguous. Despite the difficulties, we can identify a number of edge states with some confidence. 

We end the Appendix by making several observations. First, at this relatively small system size, bulk excited 
states can have energies as low as those of edge states. In fact, the lowest energy eigenstates for AM = 1-3 are 
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TABLE VIII: Overlap matrix of A = 0.5 (row) and 0.0 (column) for N = 12, N orb = 26, M = 129, and d = 0.6. The underlined 
numbers are the overlap between an edge state for A = 0.5 and the (likely) corresponding edge state for the pure Coulomb case. 
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TABLE IX: Overlap matrix of A = 0.5 (row) and 0.0 (column) for N = 12, N orb = 26, M = 130, and d = 0.6. The underlined 
numbers are the overlap between an edge state for A = 0.5 and the (likely) corresponding edge state for the pure Coulomb case. 



bulk states, with small but finite overlaps (~ 0.2) with the corresponding edge states (see Tables IVTl IVIII and Willi . 
This suggests in the pure Coulomb case fermionic edge states mix with bulk states, which is consistent with the 
fact that fermionic edge states extrapolated from finite-A neutral velocities are expected at energies in between the 
corresponding lowest two levels for AM = 2 and 3 (see Fig. [5]). A recent density-matrix renormalization group 
(DMRG) calculation 40 estimate the excitation gap to be about 0.03 e 2 /eZs, thus we expect these bulk states will float 
up in the thermodynamic limit. Second, the low-lying fermionic edge excitations do exist for A = at small excitation 
energies. However, we cannot easily decompose these states into Majorana fermionic levels with linear dispersion 
relation as we have done for A = 0.5 and 0.1. The difficulty is due to mixing of the fermionic edge excitations with 
bulk states. Third, there is significant redistribution in the weight of the lowest two edge excitations for AAf = 4 as 
A decreases, as indicated by the overlaps of the two states for A = 0.1 and 0.0 with those for A = 0.5 (see Tables IVl 
and [Dp). 
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